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Abstract 

The energies of 21 and 22 bound states of S-, P-, D-, and F-symmetry of the three-electron 
systems Li atom and Be+ ion, respectively, are calculated using the same quantum mechanical 
methods for all states, namely, the Configuration Interaction (CI) with Slater orbitals and L-S 
eigenfunctions, and the Hylleraas Configuration-Interaction (Hy-CI) methods. In this paper we 
discuss the construction and selection of the configurations in the wave functions, optimization of 
the orbital exponents and further computational techniques in order to calculate accurately the 
energy of the states using compact wave functions. The obtained accuracy for the energy using 
the Hy-CI method is about 1 x 10~^ a.u. in the case of low lying states and 1 x 10"'* a.u. for high 
excited states, with respect the nonrelativistic values. In the case of the CI method we obtain an 
overall of 1 x 10^^ a.u. accuracy. We present also a few benchmark calculations and new values 
for very high excited states. 
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I. INTRODUCTION 



Nowadays Li atom has become, like the He atomyears before, a system to test Quantum 
Chemistry and High Precision Atomic Physics 

llllal- The non-relativistic wave functions 
of three-electron atoms and ions are of great interest in applications related to highly ac- 
curate evaluations of the lowest-order relativistic and QED corrections. At this moment 
we do not have any closed procedure which can be used to construct Dirac-type (or man- 
ifestly covariant) wave functions for two- and three-electron systems. Therefore, in actual 
applications such few-electron wave functions are approximated by the solutions of the non- 
relativistic Schodinger equation(s). All corrections are evaluated with the use of the regular 
Rayleigh-Schodinger perturbation theory. It explains why the non-relativistic wave func- 
tions of three-electrons atoms and ions are of paramount importance in applications. On 
the other hand, the accuracy of modern laser-based atomic experiments allows one to deter- 
mine many transition lines (or transition energies) in three-electron atoms and ions to the 
accuracy which could be expected even twenty years ago. To match these experimental re- 
sults we need to increase (and very substantially) the accuracy of our current non-relativistic 
three-electron wave functions. 

In the last few years the low lying states of the Li atom have been calculated to the 
accuracy from a nanohartree to beyond a picohartree (1 x 10^^ — 10^^^ a.u.) 443]. These 
calculations which serve to test the Quantum Mechanical methods, are performed once, and 
their requested computational times or the length of the wave function expansions are of no 
importance. They are benchmark calculations. Usually the corresponding wave functions 
are extremely long, in the order of several ten thousands of configurations. This fact makes 
it impossible to employ these wave functions in production calculations or the calculation 
of properties, and this despite of the development of the computer technology. Therefore 
there is a need of accurate methods with accuracies in the order of a cm~^ (spectroscopical 
accuracy). For this reason, the research on compact wave functions is still in continuous 



development |8l4l2|. 



For the calculation of properties it would be desirable to have at hand all energies and 
matrix elements for ground and all excited states, and in addition, these states should be 
calculated with approximately the same accuracy. Moreover numerous excited states of all 
symmetry types (S, P, D, F, . . .) are usually needed, i.e. for the calculation of the probability 
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of ionization. Not less importantly, the computational time should be acceptable. Here the 
example of the Li atom can serve to test the methods and techniques which are developed 
for the calculations of properties, like excitation energies, transition probabilities, ioniza- 
tion energies, analysis of optical spectra, e nerg y levels in confinement conditions, nuclear 
reactions and /3^-decays, etc, see, e.g., Ref. [l3| 

In this work we employ the Hylleraas- Configuration Interaction method (Hy-CI) and 
the Configuration Interaction (CI) method with Slater orbitals and L-S eigenf unctions to 
calculate all states of the Li atom and Be"*" ion which lay below their respective energy 
limits of electronic ionization. In the next sections we will discuss the technique to select 
the energetically important configurations, and the optimization procedure of the orbital 
exponents in order to calculate accurately compact wave function expansions. Using this 
method we have obtained several benchmark energies. 



II. THE HY-CI AND CI METHODS 



The Hy-CI method was proposed by Sims and Hagstrom |14|-|lq|. The advantage of the 
Hy-CI method with respect to the other Hylleraas-type methods is that only up to one 
interelectronic coordinate rjj per configuration is introduced into the wave function and 
therefore, the method can in principle be applied to any atom. Calculations using Hy-CI 
wave functions on the atoms from He to B and the H2 molecule are reported in Refs. 
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22| . The CI wave function with Slater orbitals and L-S eigenfunctions can be considered as 
a basic part of the Hy-CI wave function. Therefore in this work we start our calculations 
with the CI wave functions. In this respect, we follow the same method as Weiss and Bunge 
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24| and use relatively small basis sets. Recent extensive CI calculations with Slater 
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26|. Both 



orbitals on Be and B atoms which are more accurate can be found in Refs. 
Hy-CI and CI wave functions can be summarized in the following expression: 

N 

^ = J2 Cp'^p^ % = 0{L'')Acl)pX (1) 

p=i 

The Hy-CI and CI wave functions are linear combinations of symmetry adapted configu- 
rations $p and the coefficients Cp are determined variationally. In this work the symmetry 
adapted configurations are constructed 'a priori' so that they are eigenfunctions of the an- 



3 



gular momentum operator LP'. Another possibility would be the posterior projection of the 
configurations over the proper spatial space, as signalized in Eq. (1) by the projection oper- 
ator 0{LP). a is teh antisymmetrization operator, and the spin eigenfunction x- In the case 
of the Li atom it is sufficient to use only one spin-function (formally a linear combination of 
the two possible spin eigenfunctions would be necessary): 



(2) 



since the energetic contribution of the second spin eigenfunction has been proven to be small 



fin the order of 1 x 10" 



i.u. 0,Q)- 



Moreover the Slater determinants produced by the 



second spin eigenfunction {2aa(3 — f3aa — af3a) (due to the antisymmetry) are repeated 
when considering already the first spin eigenfunction {af3a — (3aa). The spatial part of the 
basis functions are Hartree products of Slater orbitals: 



(3) 



k=l 



The powers z/ = 0, 1 stay for CI and Hy-CI wave functions, respectively. Powers u > 1 are 
effectively reduced to z/ = 0, 1 since all even and odd powers of be expressed as a 

product of Vij times a polynomial in r^, rj and angular functions. 

The basis functions 0p, are products of s-, p-, d- and /-Slater orbitals. Higher angular 
orbitals are in practice not needed. They are needed to obtain an accuracy of about a 



nanohartree (1 x 10 ^ a.u.) or higher (see Ref. j4|). We use unnormalized complex Slater 
orbitals for which the exponents are adjustable parameters, defined as: 



(4) 



The spherical harmonics with Condon and Shortley phase 

'2l + l{l~myy^ 
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p. 52] are given by: 
P™(cos^)e*™^, (5) 



47r (/ + m)!_ 

where P[^{cos6) are the associated Legendre functions. The spherical harmonics and asso 
ciated Legendre functions used along this work are written explicitly in {29, p. 14] 
The Hamiltonian in Hylleraas coordinates may be written 
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2 ^ TiTij dridvij 2 ^ f--^ VijVik dnjdrik 



— sin {^i - ^j) ■ (6) 

r^rij sm fc^j oipidrij 



The angular momentum operator can be extracted from Eq. (6): 



V^f2_ lyi^ 1 ^ ly^ cot^, a 

1=1 ' ' 1=1 * I. 1 1 J 



1=1 



and its eigenvalue equation used: 



= + 1)04, (8) 



with /j the angular quantum number of the orbital (pi. In the case of Hy-CI wave functions 
the term d'^ / {drijdrik) containing derivatives with respect to two rjj vanishes. 

The kinetic energy operator has been separated into several radial and angular parts. 
This operator has the advantage that for the case of three-electron kinetic integrals the 
expansion of rij into r< and r> is avoided, and therefore no three-electron auxiliary integrals 



W integrals are needed, see Ref . 32] . This fact saves not only calculations but also memory 
space. Only the easily computed two-electron auxiliary integrals V{n,m; a, (3) are needed. 
From the variational principle one obtains the matrix eigenvalue problem: 



(H - ES)C = 0, (9) 



where the matrix elements are: 



Hki = J ^kH^idr, Ski = j ^Mt. (10) 
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The occurring integrals in Hy-CI calculations of three-electron systems are two-, and three- 
electron integrals. The two-electron integrals are of the types: 



(^2), (^12), 
('^i2)(r34), (ri2) 



1 

1 

^34 



(11) 



The notation (ri2) represents the integral where the left and right hand side orbitals of 
electrons 1,2 are involved (0(ri)0(r2)ri20(ri)0(r2)). In this work the two-electron integrals 



have been calculated with the algorithms described in Ref. 
are of the following types: 



33| . The three-electron integrals 



n2ri3 



runs 



ri2 
ri3 



\ ^"23 / 



(12) 



the first three cases are evaluated by direct integration over both, one r^j and one electron, 
and thus reduced to a linear combination of two-electron integrals 3^. Only the so-called 
triangle integral (ri2ri3/r23) has not been developed by us and we use a very efficient sub- 
routine from Sims and Hagstrom [35|. Finally the two- and three-electron kinetic energy 
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33|. 



integrals have been evaluated using the Hamiltonian of Eq. (6) 

The method of integration consists in the direct integration over the interelectronic co- 
ordinates Tij and simultaneously over the coordinates of one electron. The integration of 
a three-electron integral leads then to a limited linear combination of two-electron inte- 
grals. These should be calculated very accurately in terms of two-electron auxiliary integrals 
V{rn^ n; a, /3), defined: 



l/(m, n; a, /5) = / ri'"e 
'0 



r2 



(13) 



The two-electron auxiliary integrals with positive indices m, n are evaluated efficiently in 



our computer program in terms of one-electron auxiliary integrals A{n,a) 36|. Altogether 



only two-electron integrals, as in the case of the CI method, and triangle integrals have to 
be computed. This fact will be extremely helpful when extending the application of Hy-CI 
method to larger systems. In our code, approximately the same computer memory is needed 
for CI and Hy-CI calculations. Note, that in the Hy-CI method for any electron number 
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(starting with four electrons systems) no higher order integrals than four-electron ones will 
occur. 

For our calculations in this study we have written a three-electron Hy-CI computer pro- 
gram for three-electron systems in Fortran 90. The calculations have been done with the use 
of the quadruple precision. The program has been thoroughly checked comparing results o " 
our numerical calculations with analogous results by Sims and Hagstrom 14 1 and King jo] 
for the lithium atom. Note that in such calculations we have obtained complete agreement. 

III. CALCULATIONS 

A. Construction and selection of the symmetry adapted configurations 

In this work we have performed a systematical selection of the CI and Hy-CI configu- 
rations according to their energy contribution. The ground state configuration of the Li 
atom and Be"'' ion is sss (i.e. s(l)s(2)s(3)) and the configurations for the S-symmetry states 
(L=0) that have been selected in energetic order are spp, pps, sdd, dds, sf f and ffs. The 
energetically important configurations for L = 1,2,3 are listed in Table I. The quantum 
number M=0 has been chosen because usually for this case an smaller number of Slater 
determinants is required. More types of configurations can be constructed for L=0. For 
instance, configurations like psp if the exponents ai ^ ai could be considered. In this 
work the orbital exponents in the K-shell are kept equal, see Tables III and IV, therefore 
configurations like psp are equivalent to the configurations spp. 

Other possible configurations with L=0, M=0 like ppp do exist but were discarded after 
being tested in actual calculations, when their overall energy contribution was less than 
1 X 10^^ a.u. Table II shows how the configurations used in this work have been constructed 
with S-, p-, d- and /-Slater orbitals. In addition, there are more possible 'degenerated L- 
eigenfunctions' solutions with a larger number of Slater determinants (they are degenerated 
with respect to the quantum numbers L and M but with possible different energy contribu- 
tion, i.e. non-degenerate with respect to the energy) {2^. Although the inclusion of various 
degenerated configurations has been proved to improve the energy of the state, this con- 
tribution is very small. This is important for very accurate CI calculations, like the ones 



of Bunge 
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26| . In our work, we have concentrated in the energetically most important 
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CI configurations in order to use them as a basis for Hy-CI configurations (configurations 
multiplied by an interelectronic distance rij). 

After selecting the types of configurations, we have constructed complete blocks of these 
configurations for a given basis set. For instance for the basis n=4 (i.e. [4s3p2dlf] or 
[ls,2s,3s,4s,2p,3p,3d,4f]) in the sss block we have considered the configurations: lsls2s, 
Is2s2s, 2s2s2s, IslsSs, Is2s3s, 2s2s3s, . . . , 4s4s4s. Note that the configuration Islsls has 
no physical meaning but posses a large energy contribution. Altogether our CI calculations 
can be considered 'selected' with respect to the type of configuration, and 'fuU-CI' with 
respect to the orbitals basis set. 

Another important aspect in CI and Hy-CI calculations is the symmetry adaptation of the 
configurations. The configurations arc constructed 'a priori' so that they are eigenfunctions 
of the angular momentum operator . In the sums of Table II the configurations are formed 
by Slater determinants. The determinants are pairwise symmetric (i.e. and in 

the spp configuration) and lead to the same values of the electronic integrals. Therefore it is 
possible and desirable to consider only one of the determinants and to figure out the result 
by the other. Concretely, the solutions of the eigenvalue problem, using reduced matrix 
elements 1x1 (where the integrals are added, configuration + or explicit 

matrix elements 2 x 2 of the Slater determinants, are the same. The symmetry adaptation is 
computationally favorable, since the number of Slater determinants in the input is smaller 
and repeated computation of equal integrals is avoided. By observation of table II, the 
procedure can be used for all constructed configurations. 

The Hy-CI configuration blocks have been constructed by including (1) the corresponding 
CI block; (2) the CI block multiplied by the interelectronic coordinate ri2; (3) the CI block 
multiplied by ria; and (4) the CI block multiplied by r23. Here one has to take into account 
the existing symmetries between equivalent configurations. This fact can produce linear 
dependences which make that the calculation break down (linear dependent equations in 
the eigenvalue problem). For example, 2s2s3s x ri3 is equal/equivalent to 2s2s3s x r23- 

In general, energetically important Hy-CI configurations must not be the same than the 
corresponding CI ones, but usually this is the case. Therefore wc have constructed the Hy- 
CI blocks of configurations from the selected CI ones. The number of configurations grows 
very fast when adding the three rij factors. Therefore we have filtered the configurations 
within a block one by one with the energy criterion l-Bj-i — E^] > 10~'''a.u. This means that 
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all configurations have been checked and the ones whose energy contribution is small have 
been discarded. These configurations are kept and could be included in more accurate wave 
functions. Subsequently, configurations with higher n were added and filtered. The selection 
process has been carried out for every state. Therefore the calculated wave functions for the 
Li atom and Be+ ion in ground and excited states include different configurations and have 
different lengths. 



B. Optimization of the orbital exponents 

The orbital exponents have been optimized for each atomic state of the Li atom and Be"*" 
ion. A set of two exponents has been used, one for the K-shell and the other for the odd- 
electron of the L-shell. It was kept equal for all configurations. This technique accelerates 
the computations obtaining sufficiently accurate results for the purposes of the calculation 
of properties. It is clear that for highly accurate energies beyond microhartree-accuracy, 
more fiexibility in the exponents is needed, as shown in recent calculations on the lithium 



atom with extensive optimization 
with carefully chosen exponents |4 




E^lllor in calculations with very large wave functions 



37j . The virial factor: 



X = -^ (14) 
X 2(T) ^ ' 

is used to check the quality of the wave function and guide the numerical optimization of 
the exponents in the trial wave functions. 

The optimization of two exponents at the same time, in the case of Li for all configura- 
tions, has the advantage that (being a global optimization) it is very fast, in contrast with the 
partial optimization of configurations one by one, which may take very long computational 
times. 

The optimization of the orbital exponents has been carried out with the help of a parabolic 
procedure. The orbital exponents are varied in by a step size. Three energy values are cal- 
culated and fitted to a parabola. The minimum of the parabola is calculated. Subsequently, 
this value is kept fixed and the same is done for the next exponent. The procedure consists 
of cycles and macrocycles. In the macrocycles, the step size is decreased by a given factor 
and the cycles of exponents are repeated. At every step the virial factor is calculated. The 
optimization is repeated until the energy no longer improves, and the best virial and energy 
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values agree. The optimization program is completely automatic and the exponents can be 
optimized for every state and nuclear charge. 

In Tables III and IV the optimized exponents of the CI and Hy-CI wave functions are 
given. Note that the CI exponents are in general larger than the Hy-CI ones. This is in part 
because for the Hy-CI wave functions we are using smaller basis set. For low lying states we 
have optimized a set of about 400 Hy-CI functions. For higher states (> 4) we have taken 
the CI exponents for the Hy-CI wave functions. 

IV. RESULTS 

We have calculated the S-, P-, D-, and F-symmetry states of Li atom with the CI method 
using the symmetry adapted configurations of Table I, starting with a basis set n = 4. 
The orbital exponents were optimized until the same energy minimum was obtained in two 
successive optimizations. Then these exponents were used in the next CI calculation with the 
basis n = 5, and optimized again, and so on until the basis n = 7. The optimized exponents 
of the basis n = 7 are given in Table HI. For the Be"*" ion the same procedure was repeated 
with the nuclear charge Z = 4. The excited states were determined by optimization of 
the orbital exponents for the second, third, . . . eigenvalue. Note that in strictly variational 
methods, the successive excited states are the roots of the eigenvalue problem. We have 
obtained energies that are about 1 millihartree accurate (1 x 10~^ a.u.) with respect to the 
nonrelativistic values of the literature. 

The CI calculations were carried out using double precision arithmetic in the computer 
program (about 15 decimal digits accuracy in our computer). We have calculated the energy 
of 6 S-states, 6 P-states, 5 D-states and 4 F-states. Some of these states are reported here for 
first time. Due to the large number of excited states we have chosen the notation IS for the 
ground state, 2S for the first excited and so on. The hydrogenic notation would be rather 
confusing in this context. The total energies of all 21 states of the Li atom determined in 
this study are below than the total energy of the ground state of the ""Li"*" ion, i.e. Etr ~ 
-7.279 913 412 669 305 964 914 69(5) a.u. ^:]. This total energy of the ground state in the 
two-electron "^Li"*" ion is the natural threshold energy for an arbitrary bound state in the 
three-electron Li atom. 

For the Be^ ion we have determined the total energies of twenty two bound states. 
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including seven S-states, six P-states, five D-states and four F-states. Tlie results of our 
calculations can be found in Tables V and VI. In all these calculations we have applied CI 
method. The computed energies are lower than the corresponding ionization energy of Be+ 
ion [isi. The accuracy of the calculations is also of about 1 millihartree. Note that the 
F-states calculated with the CI method are reported here for first time. 

In the Hy-CI calculations we have employed the same blocks of configurations as in the 
CI calculations, see Table I, and added blocks of these configurations multiplied by one 
interelectronic coordinate at a time, i.e. CIx(l + ri2 + ri3 + r23). For details on the selection 
of configurations, see the Calculations section. We have performed Hy-CI calculations up 
to the basis n = 6. 

It is important to mention, that in Hy-CI calculations, usually it is not necessary to 
use a basis as large as in the CI calculations. The reason is the main advantage of the 
Hy-CI method: the wave function expansion converges faster to the exact solution. This 
is due to the explicit inclusion of the interelectronic coordinate in the wave function. In 
contrast, in the CI method, the interelectronic coordinate is not explicitly considered, and 
its effect is replaced by the use of high angular momentum orbitals. Therefore in Hy-CI 
calculations high angular momentum orbitals (/ > 3) are not needed to achieve an accuracy 
of microhartree, which is the purpose of this paper. 

In this work we have not calculated F-symmetry states with the Hy-CI method. This is 
of course possible. In our computer program the electronic integrals are defined for every / 
quantum number, but the kinetic energy integrals are currently restricted to / > 2, see Ref. 
32| . The theoretical and computational implementation is somewhat cumbersome, and will 
be reported elsewhere. 

We have optimized the Hy-CI orbital exponents using a basis set n = 4 of about 400 
configurations of all types considered. Subsequently, the exponents were kept fixed for 
calculations with n = 5 and n = 6 basis sets. The orbital exponents for very high excited 
states were difficult to calculate with small basis sets (higher n needed, i.e. for a 5S states 
orbitals with n = 5, 6 are needed). Therefore the optimized CI exponents were used. The 
orbital exponents of the Hy-CI calculations are shown in Tables HI and IV. 

Very similar to the CI method, the orbital exponents were optimized for every state. In 
addition, the selection of the contributing configurations of the blocks were done separately 
for every state. Therefore, the length of the wave function and the configurations included 
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differ from state to state. This is natural, since we need different configurations to describe 
different excited states. 

The nonrelativistic energies of the first 5 S-symmetry states and the IP and ID excited 
states of Li atom and Be"*" ion are known to high accuracy, whereas the other states have 
been determined less accurately (shorter wave functions), see Tables VII and VIII. 

In the Hy-CI calculations on the Li atom, see Table VII, we have obtained an accuracy of 
1-11 microhartree (1 x 10~^) for the ground state and first S-, P-, and D-excited states. For 
higher excited states the accuracy is about 50 microhartrees. For very high excited states 
it is about 1 millihartree, as far as values for these states were known. We have determined 
the A^D -7.293 405 644 a.u. and -7.289 731 555 states for first time, and got the best 
energy to date (benchmarks) for the 3^P states -7.311 776 235 a.u. and -7.311 166 857 
a.u. 

For the Be"*" ion the accuracy is still better, probably due to the larger nuclear charge, 
see Table VIII. For the ground state and lower lying states of every symmetry we have 
obtained an accuracy of few microhartrees, using less than 1000 configurations, whereas the 
best calculations in the literature needed ten times more (up to 35 000) configurations. The 
benchmarks obtained for Be"*" ion are: the 3^P state with -13.783 519 845 a.u., and the herein 
newly calculated states: -13.699 131 127 a.u., -13.699 131 127 a.u., -13.709 
474 219 a.u. and 5^D -13.709 474 219 a.u. The dissociation threshold for the three-electron 

n 

°°Be+ ion is ^ -13.655 566 238 423 586 702 079 5(10) a.u. [38]. This value coincides with 
the total energy of the ground state of the Be^"^ ion with the infinitely heavy nucleus. 

We have obtained this accuracy employing only 0.3 % of the configurations employed in 
the most highly accurate calculations reported. Finally, all calculated states are ordered by 
their energy and presented in Tables XIX and X. All these states lay below the energy of 
ionization. Possibly, among these states, G-symmetry states with high energy are missing. 
In future works we plan to consider these states as well. 

V. CONCLUSIONS AND PERSPECTIVES 

We have calculated the total energies of twenty one and twenty two bound states in the 
Li atom and Be "*" ion, respectively. The variational wave functions of the S-, P-, D- and F- 
bound states in these three-electron atomic systems were constracted with the use of the CI 
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and Hy-CI methods. The procedure consisted in the appropriate selection of configurations 
and optimization of one set of orbital exponents for every state. The total energies of the 
low-lying states are microhartree accurate, while for excited states the accuracy is about 
1 X 10^^ — 10~^ a.u. We have obtained several benchmarks and reported for first time the 
energy of highly excited states. These wave functions are convenient for the calculation 
of properties. Concretely, these wave functions have been used for the calculation of the 
transition probabilities during nuclear /3-decay, where wave functions of very good quality 



are necessary to describe the atomic effects during nuclear reactions |45l-l47l|. We plan to 
use this systematic method of calculation to determine ground and excited states of the 
following atoms and isoelectronic ions in the periodic table, such as Be, B and C. 
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TABLE I: Symmetry Adapted Configurations employed in the calculations of the S-, P- D- and 
F-symmetry atomic states of Li atom and Bc+ ion using the CI and Hy-CI methods. The Hy-CI 
configurations are obtained from the CI ones by multiplying them by the factor R = {1 + ri2 + 
ri3 + 7-23}- 

Atomic State Symmetry adapted configurations 
Configurations 

ls2 2s "^S sss, spp, pps, sdd, dds, sff, ffs 

ls2 2p ssp, sps, ppp, ddp, pdd, spd, pds, sdp 

ls2 3d ssd, sds, spp, pps, ppd, ddd, spf, pfs, sfp 

ls2 4f ssf, sfs, ppf, ddf, flf 
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TABLE II: List of the Symmetry Adapted Configurations with quantum numbers L=0, 1,2,3 and 
Mz=0 employed in the CI and Hy-CI calculations of the Li atom and Be"*" ion. The notation sss 
stands for s(l)s(2)s(3). 



L Confs. Construction 

.sss sss 

spp spopo - spip-i - sp-ipi 

pps popos - piP-is - p-ipis 

sdd sdodo — sdid-i — sd-idi + sd2d-2 + sd-2d2 

dds dodos — did-is — d-idis + d2d-2S + d-2d2S 

aff sfofo - s/i/_i - s/-i/i + S/2/-2 + S/-2/2 - S/3/-3 - S/-3/3 

ffs fofos - /1/-1S - /-i/is + /2/-2S + /-2/2S - /3/-3S - f-sfss 

1 ssp sspo 
1 sps spos 

1 ppp PoPoPo - PiP-iPo - P-lPlPO 

1 ddp dodopo - did-ipo - d-idipo + ^2^-2^0 + d-2d2Po 

1 pdd pododo - podid-i - pod-idi + pod2d-2 + pod-2d2 

1 spd spodo — spid-i — sp-idi 

1 pds podos — pid-is — p-idis 

1 sdp sdopo — sdip-i — sd-ipi 

2 ssd ssdo 
2 sds sdos 

2 spp spoPo + spip-i + sp-ipi 

2 pps popos + pip-is + p-ipis 

2 ppd PoPodo - PiP-ido - p-ipido 

2 ddd dododo — did-ido — d-idido + ^2^-2(^0 + d-2d2do 

2 spf spofo - spif-i - sp-ifi 

2 pfs pofos - pif-is - p-ifis 

2 sfp sfopo - s/ip-i - sf-ipi 

3 ssf ssfo 
3 sfs sfos 

3 ppf PoPofo - PiP-ifo - P-mfo 

3 ddf dodofo- did-ifo- d-idifo + d2d-2fo + d-2d2fo 

3 /// /o/o/o — /1/-1/0 — /-i/i/o + /2/-2/0 + /-2/2/0 — /-3/3/0 — /-3/3/0 
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TABLE III: Exponents used in the CI and Hy-CI calculations of the Li atom. 



State 


aci = /3ci 


7CI 


aHy-CI = /^Hy-CI 


7Hy-CI 




4.644060 


1.107868 


2.994250 


0.839625 




4.698079 


0.561144 


3.550050 


0.438800 




4.605560 


0.359164 


3.241342 


0.304425 


425 


4.639431 


0.258794 


3.906990 


0.235023 


525 


4.602371 


0.191464 


4.602371 


0.191464 


625 


4.696442 


0.154444 






l2p 


4.451592 


0.827973 


3.809592 


0.762425 


22p 


4.507292 


0.504441 


3.520217 


0.370508 


32p 


4.486842 


0.320982 


3.523842 


0.255675 


42p 


4.577346 


0.230959 


3.844842 


0.204309 


52p 


4.513159 


0.169635 


3.426592 


0.170633 




4.581002 


0.144740 






l^i? 


4.512037 


0.459812 


3.359717 


0.347508 


22i? 


4.483217 


0.247508 


3.496259 


0.253341 


32£i 


4.483288 


0.200160 


3.846369 


0.337295 


4?D 


4.539008 


0.166704 


4.218008 


0.166704 


52i? 


4.492642 


0.143125 


4.492642 


0.143125 


l^F 


4.630029 


0.299372 






22f 


4.645336 


0.186989 






32F 


4.669615 


0.174730 






42F 


4.763930 


0.141984 
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TABLE IV: Exponents used in the CI and Hy-CI calculations of the Be"*" ion. 



State 


aci = /3ci 


7CI 


aHy-CI = pHy-Cl 


7Hy-CI 




6.345407 


1.950188 


4.173235 


1.406372 


2^S 


6.307744 


1.019290 


4.838717 


0.876476 




6.280890 


0.665840 


4.766466 


0.585863 




6.327154 


0.479516 


5.331658 


0.442445 




6.304912 


0.364545 


6.304912 


0.364545 




6.341861 


0.304815 


6.341861 


0.304814 




6.685487 


0.167770 






12P 


6.141069 


1.760345 


4.746625 


1.321000 


22p 


6.148030 


0.964810 


4.837058 


0.712777 


32p 


6.158844 


0.631221 


4.800125 


0.516000 


42p 


6.168223 


0.465315 


5.713473 


0.414065 


52p 


6.189834 


0.344862 


3.951567 


0.339598 


62p 


6.222913 


0.290957 


4.502617 


0.340623 


127? 


6.132817 


0.867614 


4.804284 


0.670348 


22i? 


6.159980 


0.591028 


4.750784 


0.588652 


32i? 


6.134907 


0.434209 


5.720102 


0.381441 


42i? 


6.484197 


0.337802 


6.164669 


0.296639 


52i? 


6.157877 


0.286801 


6.164653 


0.266044 


12F 


6.440871 


0.589239 






22f 


6.416382 


0.372527 






32F 


6.352375 


0.344411 






42p 


6.419948 


0.285310 
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TABLE V: Convergence of Full-CI (L-S) calculations on the ground and excited states of the Li atom with respect to the basis set. 
Basis sets constructed with s-, p-, d- and f-Slater orbitals. N is the number of configurations. The basis n = 4 stays for [4:s3p2dlf] or 
[Is2s3s4s2p3p4p3d4d4/]. Energy in a.u. and energy differences in mh. 



State 


N 


n=5 


N 


n=6 


N 


11=7 


N 


Ref. Energy 


Ref. 


diff. 




311 


-7.476 304 44 


596 


-7.476 816 76 


991 


-7.477 191 60 


34020 


-7.478 060 323 910 146 894 


\7] 


0.87 




311 


-7.352 351 26 


596 


-7.352 980 75 


991 


-7.353 248 72 


9576 


-7.354 098 421 444 364 045 


\7] 


0.85 




311 


-7.316 684 51 


596 


-7.317 410 04 


991 


-7.317 678 91 


34020 


-7.318 530 845 998 906 901 


17] 


0.85 


425 


311 


-7.301 741 72 


596 


-7.302 341 81 


991 


-7.302 681 63 


34020 


-7.303 551 579 226 734 650 


17] 


0.87 


525 






596 


-7.294 675 81 


991 


-7.294 935 22 


17072 


-7.295 859 509 943 


\£ 


0.92 


62s 










991 


-7.289 596 50 










12P 


440 


-7.407 946 09 


849 


-7.408 436 70 


1430 


-7.408 618 73 


32200 


-7.410 156 532 652 370 


[7] 


1.54 


22p 


440 


-7.334 958 50 


849 


-7.335 436 16 


1430 


-7.335 658 20 


3012 


-7.337 151 70 


[39] 


1.50 


32p 


440 


-7.309 724 51 


849 


-7.310 199 70 


1430 


-7.310 382 99 


120 


-7.311 736 


[40] 


1.34 


42p 


440 


-7.298 193 09 


849 


-7.298 615 12 


1430 


-7.298 802 13 


1021 


-7.300 287 5 


m 


1.50 


52p 






849 


-7.292 379 98 


1430 


-7.292 544 80 




-7.294 020 


[42] 


5.27 


62p 










1430 


-7.288 748 91 










l2_D 


384 


-7.333 500 02 


646 


-7.333 934 72 


1056 


-7.334 100 04 


32760 


-7.335 523 543 524 685 


[U 


1.42 


2^D 


384 


-7.309 181 71 


646 


-7.309 598 04 


1056 


-7.309 761 50 


841 


-7.311 190 1 


[43] 


1.43 


3^D 


384 


-7.297 924 73 


646 


-7.298 339 69 


1056 


-7.298 501 90 




-7.299 928 


[9] 


1.43 


A^D 






646 


-7.292 225 44 


1056 


-7.292 387 16 










5^D 










1056 


-7.288 700 67 










l^F 


120 


-7.309 517 06 


286 


-7.310 288 32 


532 


-7.310 609 58 




-7.311 168 7 


\9] 


0.56 


22f 


120 


-7.298 267 04 


286 


-7.298 989 11 


532 


-7.299 340 26 




-7.299 917 1 


[9] 


0.58 


32f 






286 


-7.292 768 68 


532 


-7.293 211 21 










42F 










532 


-7.289 401 46 











TABLE VI: Convergence of Full-CI (L-S) calculations on the ground and excited states of the Be"*" ion with respect to the basis set. The 
basis set used is the same as for the Li atom calculations of Table V. Energy in a.u. and energy differences in mh. 



State 


N 


11=5 


N 


n=6 


N 


11=7 


N 


Ref. Energy 


Ref. 


diff. 


1 2 c 

lb 


311 


1 A 000 ''7 ''7 A 1 a 

-14.SZ2 774 Id 


596 


-14.323 467 80 


991 


1 A 000 '"ICtO CC 

-14.323 768 55 


13944 


1/1 00 /I "Tf'O T "TC Tnn /1 / ooX 

-14.324 763 176 790 43(22) 


[6] 


0.99 


o2 c 

2 b 


311 


-13.920 770 83 


596 


1 r\m coo 

-13.921 528 75 


991 


1 noi oon no 

-13.921 830 02 


10000 


1 noo Ton o/^o a a o 

-13.922 789 268 544 2 


[37] 


0.96 


6 b 


oil 

311 


1 'ync ^00 on 

-13.79D d32 80 


596 


-13.797 443 52 


991 


1 '7n'7 Tco 00 

-13.797 753 82 


1 000 
1888 


-13.798 716 57 


rnl 


0.96 


i b 


311 


-13.742 568 73 


596 


-13.743 267 09 


991 


10T/10 J^C/I OA 

-13.743 654 84 


1091 


10 T-1/1 £?on /3 

-13.744 630 6 


{ A A^ 
[44| 


0.98 


b 






596 


1 T1 /I 010 CO 

-13.714 813 58 


991 


1 T1 C 000 1 T 

-13.715 222 17 


2058 


1 T1 OO/^ /I 

-13.716 286 24 


rnl 

[El 


1.1 


r2 o 
6 6 






596 


-13.689 752 76 


991 


1 /''f\'~7 /I 1 00 

-13.697 421 32 










72 c 
I b 










yyi 


1 i^jon ni^o oc 










1 2 D 

IF 


440 


-14.176 680 58 


849 


1/11 TT 01 ri no 

-14.177 zlU U3 


1430 


1/11 TT /ino no 

-14.177 408 92 


10000 


1/11 Tn 000 ono O/io t 

-14.179 333 293 342 7 


roTi 


1.92 


o2 D 


3o4 


10 000 cnc no 

-13.882 595 98 


/I n 

849 


1 000 1 01 
-lo.ooo 1/4 zl 


1 A on 

1430 


10 000 /IO/I c'y 

-13.883 424 67 




1000c 1 r: 

-13.885 15 


[42J 


1 TO 

1.73 


o2 D 


384 


-13.781 120 58 


849 


1 TO 1 '1 A A C 

-13.781 744 85 


1430 


1 T01 nT/1 TT 

-13.781 974 77 




10 TOO CIO 

-13.783 518 3 


f/l 1 1 

[41J 


1.54 


/|2 7-> 


384 


-13.734 716 26 


849 


-13.735 200 34 


1430 


-13.735 465 95 


1021 


-13.737 18 


r ^ oi 

[42] 


1.71 


c2 o 
r 








1 Ti n 1 on c /I 
-lo. / lU loy 04 


1 A on 
14oU 


1 T1 n ooA on 
-lo. ( lU ooU oU 




1 T1 n^; 
-lo. Iz Ud 


[42J 


1 TO 


e.2 D 
D r 










1 A on 
14oU 


1 cnc 007 nn 
-1o.D9o zz / 99 










1 2 n 


00 /I 


1 OTC ^0 1 /in 
-lo.o75 781 49 


d4d 


1 OTC Oi^n Tn 
-1o.o7d zdU 79 


1 ncc 


10 OTC A Ad TO 

-13.o7d 44d 73 


/I 1 

o41 


1 OTO nco c 
-13.o7o (Joz d 


[43| 


1 £!n 






-±o. ( ( 4:Zo yo 




-10. ( ( oyu iZ( 




1 77Q ns'^ 

-±0. ( ( a UOO 00 


04-L 


-±0. 1 OU Uc/U u 








384 


-13.733 383 90 


646 


-13.733 841 36 


1056 


-13.734 023 64 


841 


-13.735 637 


[43] 


1.61 








646 


-13.709 376 51 


1056 


-13.709 537 57 




















1056 


-13.694 803 96 












120 


-13.778 730 68 


286 


-13.779 403 25 


532 


-13.779 945 95 












120 


-13.733 732 72 


286 


-13.734 363 32 


532 


-13.734 923 61 
















286 


-13.709 787 93 


532 


-13.710 457 22 










42^ 










532 


-13.695 578 93 











TABLE VII: Calculated Hy-CI energies of the ground S-state and first S-, P-, and D-excited states of Li atom. Convergence of the calculations 
and comparison with energy values of the literature. N is the number of symmetry adapted configurations. Energy in a.u. Energy differences 
in /ih. Benchmark calculations are written in bold. 



11 


State 


N 


n=4 


N 


n=5,6 


N Ref. Ener. 


Ref. Diff. in \jh 


1 




309 


-7.478 053 222 500 


693 


-7.478 058 968 847 


34020 


-7.478 060 323 910 146 894 


\^\ 


1.41 


3 




307 


-7.354 078 274 550 


549 


-7.354 093 706 422 


34020 


-7.354 098 421 444 364 045 


[7] 


4.71 


6 


o9 n 


252 


-7.318 481 008 275 


591 


-7.318 517 758 979 


34020 


-7.318 530 845 998 906 901 


m 


13.1 


10 


a9 n 
4'^S 






687 


-7.303 496 699 422 


34020 


-7.303 551 579 226 734 650 


[7] 


60.7 


14 








491 


-7.295 739 602 768 


17072 


-7.295 859 509 943 




120.0 


o 
z 


1 2 p 


ooi 


-/.41U lo4 IZZ DDo 


DiD 


-(AW i4y Uoy Dyo 


32200 


-7.410 156 532 652 370 


m 


7.5 


4 


22p 


530 


-7.337 055 166 688 


766 


-7.337 099 254 449 


3012 


-7.337 151 70 


[39] 


52.5 


7 


32p 


466 


-7.311 724 861 462 


659 


-7.311 776 235 496 


120 


-7.311 736 


[40] 


-40.0 


12 


42p 






573 


-7.300 091 953 421 


1021 


-7.300 287 5 


m 


195.5 


16 


5^P 






1297 


-7.293 967 121 805 




-7.294 020 


[42] 


53.0 


5 


l^D 


188 


-7.335 505 134 976 


238 


-7.335 511 801 876 


32760 


-7.335 523 543 524 685 


[7] 


11.7 


8 


2^D 


176 


-7.311 192 543 434 


187 


-7.311 211 047 253 


841 


-7.311 190 1 


[43] 


-21.0 


13 




273 


-7.298 186 482 327 


372 


-7.299 416 502 739 




-7.299 928 


[9] 


51.2 


17 








271 


-7.297 145 785 698 










20 








525 


-7.289 731 554 753 











TABLE VIII: Calculated Hy-CI energies of the ground S-state and first S-, P-, and D-excited states of the Be"*" ion. Convergence of the 
calculations and comparison with energy values of the literature. N is the number of symmetry adapted configurations. Energy in a.u. 
Energy differences in ^h. Benchmark calculations are written in bold. 

II State N n=4 N n=5,6 N Ref. Ener. Ref. diff 

1 l^S* 514 -14.324 757 377 542 1028 -14.324 761 678 324 13944 -14.324 763 176 790 150 [6] 1.5 
3 2^5 502 -13.922 759 980 099 1199 -13.922 784 968 156 10000 -13.922 789 268 554 2 [37] 4.3 
6 409 -13.798 520 453 185 757 -13.798 706 849 282 1888 -13.798 716 57 [9] 9.7 

10 4^5 1020 -13.744 518 213 729 1091 -13.744 630 6 [44] 112.4 

14 525* 1066 -13.716 152 057 533 2058 -13.716 286 24 [9] 134.2 

^ 18 771 -13.699 131 126 904 



2 l^P 373 -14.179 314 874 884 -14.179 326 617 395 10000 -14.179 333 293 342 7 [37] 7.2 

4 22p 499 -13.885 035 680 428 707 -13.885 072 810 365 -13.885 15 [42] 77.0 
7 3^^ 352 -13.783 432 326 337 634 -13.783 519 845 223 1021 -13.783 518 3 [41] -1.5 

11 42p 480 -13.734 833 028 296 -13.737 18 [42] 2350.0 

16 1232 -13.711 935 267 558 -13.712 06 [42] 125.0 

20 459 -13.711 378 664 909 

5 I'^D 265 -13.878 005 890 429 338 -13.878 031 745 646 841 -13.878 052 6 [43] 21.0 
9 230 -13.779 724 787 607 298 -13.780 335 992 025 841 -13.780 696 [43] 360.0 

13 250 -13.728 658 182 159 216 -13.735 055 937 333 841 -13.735 637 [43] 581.1 

17 i'^D 239 -13.709 739 027 953 

21 192 -13.677 409 085 288 



TABLE IX: The total energies of some S-, P-, and D-states of the Li atom (in a.u.) ordered by 
their energy. The total energy of ionized Li+ ions is -7.279 913 412 669 305 964 914 69(5) a.u. [38 1 
(threshold energy). 



n 


State 


E(FCI) 


E(Hy-CI) 


Ref. Energ 


;y 






Ref. 


1 




-7.477 


191 


60 


-7.478 058 91 


-7.478 060 


323 


910 


146 894 


m 


2 


12P 


-7.408 


618 


73 


-7.410 149 07 


-7.410 156 


532 


652 


370 


m 


3 




-7.353 


248 


72 


-7.354 093 71 


-7.354 098 


421 


444 


364 045 


m 


4 


22p 


-7.335 


658 


20 


-7.337 092 97 


-7.337 151 


70 






m 


5 


l^D 


-7.334 


100 


04 


-7.335 505 13 


-7.335 523 


543 


524 


685 


m 


6 


3^S 


-7.317 


678 


91 


-7.318 517 76 


-7.318 530 


845 


998 


906 901 


m 


7 


32p 


-7.310 


382 


99 


-7.311 776 23 


-7.311 728 


51 






m 


8 


2^D 


-7.309 


761 


50 


-7.311 211 05 


-7.311 190 


1 






m 


9 


l^F 


-7.309 


517 


06 














10 


4^S 


-7.302 


681 


63 


-7.303 490 87 


-7.303 551 


579 


226 


734 650 


m 


11 




-7.299 


340 


26 














12 


42p 


-7.298 


802 


13 


-7.300 010 05 


-7.300 287 


5 






m 


13 


3^D 


-7.298 


501 


90 


-7.296 927 03 


-7.299 928 








[9] 


14 


b^S 


-7.294 


935 


22 


-7.295 735 01 


-7.295 859 


509 


943 




m 


15 


3^F 


-7.293 


211 


21 














16 


5^P 


-7.292 


544 


80 


-7.293 967 12 


-7.294 020 








[42] 


17 


4?D 


-7.292 


387 


16 


-7.292 225 21 












18 


Q'S 


-7.289 


596 


50 














19 


A^F 


-7.289 


401 


46 














20 


62p 


-7.288 


748 


91 














21 


5^D 


-7.288 


700 


67 
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TABLE X: The total energies of some S-, P-, and D-states of the Be+ ion ordered by their energy. 
The total energy of the ground state in the Be^^ ion is: -13.655 566 238 4235 867 020 795(10) a.u. 
38l | (threshold energy). 



n State E(FCI) E(Hy-CI) Ref. Energy Ref. 



1 


1^5 


-14 


323 


768 


55 


-14 


324 


761 


68 


-14 


324 


763 176 790 150 


[6] 


2 


l2p 


-14 


177 


408 


92 


-14 


179 


326 


07 


-14 


179 


333 293 342 7 


[33 


3 


2^5 


-13 


921 


830 


02 


-13 


922 


784 


97 


-13 


922 


789 268 554 2 


[33 


4 


22p 


-13 


883 


424 


67 


-13 


885 


072 


81 


-13 


885 


15 


[42j 


5 


l^D 


-13 


876 


446 


73 


-13 


878 


005 


89 


-13 


878 


052 6 


m 


6 


325 


-13 


797 


753 


82 


-13 


798 


706 


85 


-13 


798 


716 57 


m 


7 


32p 


-13 


781 


974 


77 


-13 


783 


519 


84 


-13 


783 


518 3 


m 


8 


IT 


-13 


779 


945 


95 


















9 




-13 


779 


083 


58 


-13 


780 


335 


99 


-13 


780 


696 


m 


10 


425 


-13 


743 


654 


84 


-13 


744 


589 


13 


-13 


744 


630 6 


[43j 


11 


42p 


-13 


735 


465 


95 


-13 


733 


901 


88 


-13 


737 


18 


[42j 


12 




-13 


734 


923 


61 


















13 


3^D 


-13 


734 


023 


64 


-13 


728 


692 


76 


-13 


735 


637 


m 


14 


525 


-13 


715 


222 


17 


-13 


716 


152 


06 


-13 


716 


286 24 


m 


15 


3^F 


-13 


710 


457 


22 


















16 


5^P 


-13 


710 


139 


54 


-13 


711 


935 


27 


-13 


712 


06 


m 


17 


4?D 


-13 


709 


537 


57 


-13 


728 


692 


76 










18 


6^5 


-13 


697 


421 


32 


-13 


699 


131 


13 










19 


A^F 


-13 


695 


578 


93 


















20 


6^P 


-13 


695 


227 


99 


-13 


711 


378 


66 










21 


b^D 


-13 


694 


803 


96 


-12 


631 


500 


84 










22 


725 


-13 


680 


062 


35 
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